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ABSTRACT 


Wave  breaking  is  the  only  souree  of  energy  dissipation  in  the  Battjes  and  Janssen 
(1978)  wave  transformation  model,  whieh  is  parameterized  by  a  breaking  wave 
parameter,  y .  The  Battjes  and  Janssen  (1978)  wave  transformation  model  was  calibrated 
by  Battjes  and  Stive  (1985)  and  the  calibration  was  refined  by  Morris  et  al.  (2001)  for 
waves  over  shallow  sloping  beaches.  The  objective  of  this  study  was  to  further  refine  the 
calibration  to  include  steep  beaches  for  a  range  of  wave  conditions  by  analyzing  data 
from  a  nearshore  experiment  at  Sand  City,  California.  Waves  were  measured  by  a  cross¬ 
shore  array  of  nine  pressure  sensors.  The  pressure  data  were  analyzed  for  Hrms  and 
compared  with  calculated  Hrms  by  the  model.  Results  were  largely  inconclusive,  which  is 
attributed  to  wave  reflection  from  the  steep  beach,  something  not  accounted  for  in  the 
model.  Excluding  data  collected  at  low  tides  and  allowing  the  model  to  account  for 
reflection  would  likely  reveal  a  more  interesting  outcome. 
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EXECUTIVE  SUMMARY 


The  nearshore  wave  transformation  model  by  Battjes  and  Janssen,  (1978)  is  part 
of  the  Delft-3D  nearshore  wave  and  current  model  being  considered  for  rapid  transition 
for  use  as  the  operational  model.  Wave  breaking  is  the  only  source  of  energy  dissipation 
in  the  model,  which  is  parameterized  by  a  breaking  wave  parameter,  y .  The  Battjes  and 
Janssen  (1978)  model  was  calibrated  by  Battjes  and  Stive  (1985)  and  the  calibration  was 
refined  by  Morris  et  al.  (2001)  for  waves  over  shallow  sloping  beaches.  The  objective  of 
this  paper  was  to  further  refine  the  calibration  to  include  steep  beaches  for  a  range  of 
wave  conditions  by  analyzing  data  from  a  nearshore  experiment  at  Sand  City,  California. 
Waves  were  measured  by  a  cross-shore  array  of  nine  pressure  sensors.  The  pressure  data 
were  analyzed  for  Hnns  and  compared  with  calculated  Hrms  by  the  model.  Results  were 
largely  inconclusive,  which  is  attributed  to  wave  reflection  from  the  steep  beach, 
something  not  accounted  for  in  the  model.  Excluding  data  collected  at  low  tides  and 
allowing  the  model  to  account  for  reflection  would  likely  reveal  a  more  interesting 


outcome. 
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I.  INTRODUCTION 


Prediction  of  sea  surface  elevation  (SSE)  due  to  breaking  waves  is  an  area  of 
critical  interest  in  a  variety  of  fields  including  municipal  planning  and  military 
operations.  States,  counties,  and  cities  benefit  from  a  reliable  method  of  forecasting  surf 
conditions  in  many  areas  such  as  safety  of  beach-goers  and  zoning  restrictions.  Military 
planners,  especially  Navy  and  Marine  Corps  personnel  conducting  amphibious  operations 
and  special  operations,  require  accurate  predictions  of  surf  conditions. 

The  Battjes  and  Janssen  (1978)  waveheight  model  for  surf  prediction  (hereafter 
referred  to  as  BJ78)  is  compared  with  data  acquired  on  a  steep  beach.  The  model  has 
been  successfully  tested  with  data  from  beaches  with  gentle  slopes.  (Battjes  and  Stive, 
1985)  There  are  no  known  results  of  testing  this  model  on  steep  beach  conditions.  The 
BJ78  model  is  the  basis  of  wave  prediction  in  the  surf  zone  for  the  wave,  current,  and 
morphology  model,  DELFT3D,  currently  under  consideration  for  transition  to  operational 
use  for  the  U.S.  Navy. 

The  model  was  first  calibrated  by  Battjes  and  Stive  (1985)  for  data  having 
relatively  steep  waves  and  gentle  beach  slopes.  The  model  has  a  breaking  wave 
parameter  (y),  which  was  shown  to  be  dependent  on  deep-water  wave  slope  (5^)  as 
given  by 

y  =  0.5 -I- 0.4tanh335o  (1) 

where  Sq  =  H^j Lq  in  which  is  the  offshore  wave  height,  and  Lq  is  the  offshore  wave 
length. 

Morris  et  al.  (2001)  expanded  the  calibration  to  include  waves  of  low  wave 
steepness  (swell)  and  found  an  additional  dependence  on  beach  slope  as  given  by 

y  =  0.2-l-0.32tanh^  for^g  <0.002,  (2) 

where  ^  is  the  Iribarren  number  given  by 


tan  ^ 


(3) 


1 


where  tan  p  is  the  beaeh  slope.  The  wave  and  beach  geometry  based  on  the  Iribarren 
number  for  all  of  the  comprehensive,  large-scale  nearshore  field  experiments  is 
summarized  in  Figure  1.  The  data  analyzed  here  is  from  the  RIPEX  experiment,  which 
has  not  been  previously  analyzed.  RIPEX  data  represent  large  Iribarren  numbers  for 
large  beach  slopes  and  swell  conditions. 


Iribarren  Number,  =  tan(pJ/(H^/L^)'^ 


Io9l0(g 

Eigure  1 .  Iribarren  numbers  for  conditions  used  to  calibrate  BJ78. 


The  objective  of  this  paper  is  to  calibrate  the  Battjes  and  Janssen  (1978)  model  for 
steep  beach  conditions  by  finding  the  optimal  y  values  for  model  comparison  with  data. 
The  next  section  of  this  paper  discusses  the  experiment  including  location,  time  for  which 
data  were  collected,  sensors  used  to  collect  data,  and  conditions  of  the  beach  and  surf 
Eollowing  that  is  a  description  of  the  model  and  the  methods  used  in  this  paper  to  analyze 
the  data.  The  last  two  sections  are  results  and  conclusions. 
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II.  EXPERIMENT 


The  Rip-Current  Experiment  (RIPEX)  was  conducted  from  10  April  to  22  May 
2001  on  a  steep  beach  in  Sand  City,  California.  The  beach  steepness  varied  from  1:20 
offshore  graduating  into  a  tidal  plateau  finally  rising  into  a  steep  1:5  beach  face. 
Significant  wave  heights  ranged  from  0.5  to  2m  with  periods  ranging  from  5  to  IVsec. 
(Eigure  2)  The  waves  approached  at  near  normal  incidence  to  the  local  shoreline  due  to 
severe  refraction  and  a  relatively  narrow  aperture  caused  by  the  sheltering  headlands  of 
Monterey  Bay. 


Eigure  2.  Time  series  of  H^ms  (top)  and  period  (bottom) 


Bathymetry  data  were  collected  from  the  base  of  the  dunes  out  to  a  depth  of  1 8m 
using  three  methods  to  obtain  full  coverage.  On  the  beach,  an  all  terrain  vehicle  (ATV) 
drove  over  the  dunes  with  a  Differential  Global  Positioning  System  (DGPS)  receiver. 
Eatitude,  longitude,  and  elevation  were  recorded  with  horizontal  and  vertical  accuracies 
with  6  cm  rms  error.  A  jet  ski  with  an  echo  sounder  and  DGPS  receiver  measured 
bathymetry  from  offshore  to  a  depth  of  about  Im.  Prom  the  beach  seaward  to  a  water 

3 


depth  of  approximately  Im,  a  person  with  a  DGPS  backpaek  walked  the  nearshore. 
Bathymetry  was  measured  on  nine  days  during  the  experiment. 

For  data  considered  here,  waves  were  measured  with  a  cross-shore  array  of  nine 
pressure  sensors  extending  approximately  235m  offshore  to  a  depth  of  approximately  6m. 
In  addition,  an  offshore  wave  buoy  was  located  642m  offshore  in  16m  of  water  at 
36.61955N,  121. 85921 W.  (Figure  3) 

Pressure  sensors  over  the  shoal  were  attached  to  poles  and  sunk  vertically  into  the 
sand  to  depths  of  0.5m  to  1.0m.  The  pressure  sensors  were  a  silicon-on-sapphire 
submersible  pressure  transducer  made  by  Hydracon  Co.,  model  number  2450-301.  The 
pressure  range  was  0-25psia.  The  two  furthest  offshore  pressure  sensors  were 
incorporated  within  SonTek  5MHz  ADV  Ocean  Probes. 


X(m) 

Figure  3.  Bathymetry  and  instrument  locations  of  RIPEX  experiment.  (Coordinates  are 
centered  at  36.61629093N,  121.85327237W.)  From  Renters  et  ah,  (2002) 
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The  wave  buoy  is  a  Direetional  Waverider  Mark  II,  whieh  measures  wave  height 
and  direetion.  It  is  equipped  with  a  heave-pitch-roll  sensor  Hippy-40.  Directional  data 
are  accurate  to  within  1.5°  with  an  error  of  0.5°.  Heave  resolution  is  1cm  and  accurate  to 
3%  of  the  measured  value.  Frequency  resolution  is  0.005Hz  from  0.025  to  O.lHz  and 
O.OlHz  from  0.1  to  0.59Hz.  The  high  frequency  cut-off  is  0.6Hz. 
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III.  ANALYSIS 


A.  WAVE  TRANSFORMATION  MODEL 

The  transformation  of  waves  from  offshore  to  the  beach  is  described  using  the 
energy  flux  balance  model  accounting  for  random  wave  breaking  by  Battjes  and  Janssen 
(1978).  Assuming  straight  and  parallel  depth  contours  and  stationary  wave  conditions, 
the  energy  flux  balance  is  represented  by 

dEc„  cos  6* 

- ^ - +  Z)  =  0  (4) 

dx 


where  E  is  the  mean  wave  energy,  6  is  the  angle  of  the  waves  with  respect  to  shore 
normal,  x  is  the  horizontal  coordinate  perpendicular  to  the  shoreline,  and  D  is  the  time- 
mean  dissipated  power  per  unit  area,  and  is  the  group  velocity  described  by  linear 

wave  theory: 


c 


g 


2^r, ,  a  y 

k  f  ^  sinh  2kh  )\ 


(5) 


where  h  is  depth,  k  is  wavenumber,  and  /  is  mean  wave  frequency  of  the  assumedly 
narrow  banded  spectra.  For  simplicity  it  is  further  assumed  that  the  only  source  of 
dissipation  is  wave  breaking. 

In  deep  water,  random  waves  are  well  described  as  having  a  Rayleigh  wave  height 
distribution,  which  is  a  single  parameter  distribution  in  terms  of  the  root  mean  square 
wave  height  Assuming  linear  wave  theory,  the  mean  energy  density  (£')  for  a 

Rayleigh  wave  height  distribution  is  given  by 

E^yUpgHl,.  (6) 

To  describe  wave  breaking,  it  is  assumed  that  if  a  random  series  of  waves  enters  a 
given  water  depth  (/z),  waves  larger  than  a  certain  maximum  height  (Elm)  will  break  and 
all  waves  shoreward  of  h  will  be  equal  in  height  or  smaller  than  H ^  .  These  waves  are 
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described  in  terms  of  the  clipped  Rayleigh  distribution  of  wave  heights  with  a  cutoff  at 
determined  by  h  as  given  by 


F{H)  =  1  -  exp 

V 


for  0  <  //  <  //„ 


=  1 


for  H^<H. 


(V) 


The  probability  that  at  any  given  depth  a  wave  height  is  from  a  broken  or 
breaking  wave  is  obtained  from  Equation  (7),  as 


i-g  _ 

Ina  I  J 


(8) 


where  gives  the  fraction  of  waves  that  are  broken  or  breaking  at  any  given  depth.  In 
deep  water,  ^  0  .  In  shallow  water,  Hms  increases  and,  Qi,  ^  I  implying  that  all 
waves  have  broken  or  are  breaking,  i.e.  saturation.  The  breaking  wave  criterion  (Miche, 
1954)  is  assumed 

=  0.88^  '  tanh(k/2y/0.88) .  (9) 


where  y  is  a  coefficient  that  controls  the  fraction  of  breaking  waves.  Asy  increases,  H ^ 

for  a  given  depth  increases,  decreasing  the  amount  of  dissipation  occurring  at  that 
location.  This  allows  breaking  to  begin  farther  inshore.  In  shallow  water.  Equation  (9) 
reduces  to  =  yh . 


Energy  dissipation  in  a  broken  wave(E)  )  is  assumed  to  be  similar  to  that  of  a 
linear  bore.  If  the  waves  are  periodic  with  mean  frequency  / ,  then,  for  a  single  wave, 
the  average  power  dissipated  per  unit  area  due  to  breaking  is  given  as 


D  = 


L 


fD' 


fD'  1  - 

h 


(10) 
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where  L  is  wavelength  and  e  is  phase  speed. 

To  apply  Equation  (10)  to  random  waves,  it  must  be  applied  to  broken  waves 
only,  where  the  broken  waves  are  assumed  to  have  a  height  equal  to  Hm  and  probability 
of  oeeurrenee  of  Qt,.  Sinee  most  breaking  oeeurs  where  H ^  =  h ,  H^Ai  is  dropped  from 
the  order  of  magnitude  relationship  and  Equation  (10)  beeomes 

D  =  ^QjpgHl  (11) 


where  a  is  a  eonstant  of  order  one  that  eontrols  the  level  of  energy  dissipation.  By 
eombining  Equations  (5)  and  (11),  the  rate  of  energy  dissipated  due  to  breaking  waves 
ean  be  expressed  in  terms  of  ,  depth  {h),  and  eonstants. 


The  transformation  of  H is  obtained  by  integrating  Equation  (4)  with  known 
wave  and  depth  profdes.  In  addition,  the  wave-indueed  time  average  water  level  ehange 
(77)  is  ealeulated  by  solving  the  mean  x-momentum  flux  balanee  as  given  by 


,  dri  ^ 
— —  +  pgh—  =  0 
dx  dx 


(12) 


where  the  total  depth  is  given  as  the  sum  of  the  still  water  depth,  d,  and  set-up  (-down) 


h  =  d  +  rj 


(13) 


The  eross-shore  momentum  flux  for  normally  ineident  waves  (also  ealled  radiation  stress) 
is  given  by  linear  wave  theory  as 


S..  = 


1 


2kh 


- h 

v2  smh.2kh 


E. 


(14) 


Hrms  and  peak  period  at  ADV  #2  for  a  given  hour  are  used  as  initial  eonditions. 
Setup  is  assumed  to  be  zero  at  ADV  #2.  The  model  caleulates  eonditions  every  2.5m 
starting  one  step  shoreward  of  ADV  #2.  Eor  eaeh  point’s  ealeulations,  setup  is  initially 
assumed  to  be  equal  to  the  final  setup  of  the  previous  point.  The  resulting  ealeulated 
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setup  is  then  added  to  depth  and  a  ehange  in  setup  is  calculated.  The  model  is  iteratively 
run  until  the  change  in  setup  is  less  than  0.1%  of  depth,  which  includes  bathymetry,  tides 
and  setup. 


B,  WAVE  HEIGHT  DATA 

Wave  heights  are  determined  from  pressure  sensor  data  from  April  25  to  May  21, 
2001.  Data  sampled  originally  sampled  at  8Hz  was  decimated  to  2Hz  and  separated  into 
approximately  one-hour  segments  and  checked  for  non-periodic  sampling  times  and  gaps. 
The  longest  gap  found  in  the  data  was  1 1  seconds.  Linear  interpolation  was  used  to  fill  in 
the  gaps. 

The  data  were  then  checked  for  evidence  of  vertical  movement  of  the  pressure 
sensors.  Vertical  movement  occurred  when  the  sand  about  the  mounting  pipe  liquefied 
owing  to  pipe  vibration  caused  by  pounding  of  the  breaking  waves.  Vertical  movement 
of  the  sensors  was  a  concern  especially  for  the  sensors  in  shallow  water.  To  check  this, 
pressure  sensor  data  were  examined  in  two  segments.  The  first  was  from  April  25  to 
May  5.  The  second  included  data  from  May  5  to  May  21.  In  each  segment,  a  mean 
pressure  was  calculated  over  approximately  every  50  seconds  of  data  to  fdter  out  wave 
signature.  The  graphs  of  the  resulting  tidal  variations  were  visually  examined  for  two 
types  of  shifts:  sudden  vertical  jumps  corresponding  to  an  equivalent  of  more  than  10 
percent  of  the  water  depth  and  gradual  increasing  or  decreasing  trends.  No  trends  were 
found  and  only  four  sudden  jumps  of  this  magnitude  were  found. 

To  correct  the  data  for  the  sudden  jumps,  mean  seawater  level  was  calculated 
prior  to  and  after  the  sudden  jumps  and  the  difference  added  or  subtracted  to  the  post¬ 
jump  data.  This  brought  all  data  points  from  a  given  sensor  to  the  same  reference  level. 
This  type  of  shift  occurred  in  the  three  pressure  sensors  closest  to  shore. 

Sea  surface  elevation  spectra  (SSE)  were  calculated  from  pressure  spectra  using 
linear  wave  theory  and  applying  one  of  two  transfer  functions.  The  transfer  function  for 
buried  pressure  sensors  is  given  by  Raubenheimer  et  ah,  (1998) 

H{f)  =  e’^  co^h{kh).  (15) 
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where  z  is  the  depth  (positive  downward)  the  sensor  is  buried  in  the  sediment.  For  the 
ADV  pressure  sensors  and  at  times  when  other  sensors  were  not  buried,  the  transfer 
funetion  given  by 


H(f)  = 


eos 


h{kh) 


h{k{h  +  z)) 


eos 


(16) 


was  used  to  ealeulate  SSE.  If  the  depth  of  water  over  a  pressure  sensor  was  less  than 
5em,  no  Hrms  was  ealeulated  for  that  or  any  sensor  eloser  to  the  beaeh. 


The  Hrms  was  ealeulated  from  the  SSE  spectra,  5*  (/) ,  using 


H  =2V2( 


<7 


(17) 


where  a  =  S^{f)df  . 


The  deep-water  waves  were  calculated  by  back-shoaling  the  wave  height 

measured  at  the  buoy,  H ,  to  calculate  Sq  , 


U  _  ZIT 

rmsn  rmsi, 


A' 

J 


(18) 


where  is  group  velocity  calculated  using  linear  theory  in  Equation  5. 

Measured  bathymetry  data  were  used  to  linearly  interpolate  bathymetry  for  days 
where  no  bathymetry  data  existed. 


C.  COMPARISONS 

An  optimal  value  fory  was  determined  for  each  one-hour  run  that  gave  the  least 
error  between  the  measured  and  modeled  data.  The  error  between  the  measured  Hrms  and 
each  modeled  Hrms  were  calculated  for  a  range  of  y  between  0.25  and  0.9  in  steps  of 
0.01,  where  the  error  was  calculated  as 


meas 


(19) 
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where  N  is  the  number  of  sensors  in  a  run.  If  one  of  the  sensors  in  either  set  of  data  had 
no  value,  that  differenee  was  ignored.  Onee  the  error  was  ealculated  for  eaeh  y  run,  an 
optimal  y  was  determined  by  the  smallest  error. 
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IV.  RESULTS 


Wave  breaking  was  generally  initiated  approximately  180-200m  offshore  on  the 
offshore  slope  depending  on  the  tide  elevation  and  the  wave  height.  Onee  waves  began 
breaking,  they  continued  to  break  until  they  reached  the  beach.  However,  under  high  tide 
and  small  wave  conditions,  waves  did  not  break  until  they  reached  the  beach.  (See 
Figures  4  through  7)  Data  from  pressure  sensors  1  through  3  and  6  showed  many  gaps 
due  to  exposure  at  low  tides. 


May  2,2001-1300:  Big  waves,  High  tide.  Optimized  gamma=0.74 


40  60  80  100  120  140  160  180  200  220  240 


Figure  4.  Modeled  (solid  line)  versus  measured  (dashed  line)  Hrms  with  sensor  locations 
(stars),  setup,  and  depth  (including  tide)  for  big  wave,  high  tide  conditions 
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May  3,2001-0700:  Big  waves,  Low  tide.  Optimized  gamma*0.77 


Figure  5.  Same  as  Figure  4  but  for  low  tide 


May  13,  2001-0900:  Small  waves.  High  tide.  Optimized  gamma=0  67 


40  60  80  100  120  140  160  180  200  220  240 


Distance  offshore  (m) 


Figure  6.  Same  as  Figure  4  but  for  small  waves 
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May  13,  2001-1600:  Small  waves.  Low  tide,  Optimized  gamma=0.69 


Figure  7.  Same  as  Figure  4  but  for  small  wave,  low  tide  conditions 


Time  (yearday) 


Figure  8.  Increase  in  ADV  #2  Hnns  (solid  line)  over  buoy  Hrms  (dashed  line)  on  day  136 

(May  16) 
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The  modeled  data  agreed  well  with  measured  data,  although  the  model  sometimes 
did  not  model  shoaling  well.  All  errors  assoeiated  with  optimized  gammas  were  less  than 
8  pereent  and  the  average  error  was  3  percent.  At  times  shoaling  was  over-predicted 
(Figure  5)  and  at  times  under-predicted  (Figure  9). 


May  15,  2001-0700:  Small  waves,  Low  tide.  Optimized  gamma=0.9 


40  60  80  100  120  140  160  180  200  220  240 

Distance  offshore  (m) 


Figure  9.  Modeled  (solid  line)  versus  measured  (dashed  line)  Hrms  with  sensor  locations 
(stars),  setup,  and  depth  (including  tide)  for  small  wave,  low  tide  conditions  showing 

under  forecast  shoaling 


The  optimized  values  of  y  are  plotted  versus  Sq  and  Iribarren  Number,  (Figures 

10  and  11).  While  the  error  was  small,  optimized  gamma  was  scattered.  The  calibration 
curves  as  given  by  Equations  (1)  and  (2)  are  plotted  on  each  graph  respectively  for 
comparison  purposes.  Figures  10  and  11  show  optimized  gamma  apparently  reaching  a 
maximum  limit  at  0.9.  This  is  an  artifact  of  limiting  gamma  to  less  than  0.9  which  was 
considered  greater  than  optimal  gamma  under  any  conditions. 
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Figure  10.  Optimized  gamma  versus  deep  water  wave  steepness  (so)  overlain  by  BJ85  eurve 

for  optimal  gamma  (Equation  1) 


Figure  11.  Optimized  gamma  versus  Iribarren  Number  (^)  overlain  by  Morris’  curve  for 

optimal  gamma  (Equation  2) 
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V.  DISCUSSION 


The  comparison  of  optimal  as  a  function  of  both  deep  water  wave  steepness,  so, 
and  Iribarren  Number,  show  an  unexpectedly  large  scatter.  (Figures  10  and  11)  The 
several  points  in  Figure  10  with  so  between  0.02  and  0.025  that  are  in  good  agreement 
with  the  BJ85  curve  were  compared  with  tides  and  wave  height  for  some  point  of 
commonality  but  none  was  found. 

No  significant  trends  are  seen  in  Figure  1 1 .  The  significance  of  this  is  the  same  as 
described  previously.  However,  since  Iribarren  Numbers  were  calculated  with  a  beach 
slope  of  1:5,  these  numbers  may  be  incorrect  at  low  tide  when  the  tidal  plateau  is 
exposed.  Time  constraints  prevented  re-analyzing  the  measured  data  using  only  data 
taken  at  high  tide.  The  models  inability  to  handle  shoaling  consistently  may  be  due  to 
reflection  from  the  beach.  During  some  low  tides,  the  tidal  plateau  was  exposed. 
Reflection  increases  with  beach  slope.  Reflection  from  the  steep  1:5  beach  was 
intermittent  and  occurred  only  on  higher  tides.  Since  the  model  does  not  account  for 
reflection,  this  is  a  source  of  error.  It  was  usually  observed  that  the  reflected  and  incident 
waves  interacted  to  cause  wave  breaking  when  they  collided  with  each  other.  Reflected 
waves  would  be  expected  to  cause  standing  wave  patterns  that  would  modify  how  waves 
break.  An  example  of  small  waves  at  high  tide  is  shown  in  Figure  6  when  strong 
reflection  would  be  expected  from  the  beach  face,  showing  large  differences  between 
measurements  and  the  model 

ADV  #2  gave  Hms  larger  than  that  of  the  seaward  buoy  after  May  16.  (Figure  8) 
ADV  #2’s  increase  in  Hrms  over  that  of  the  buoy  was  puzzling.  A  visual  check  of  the 
pressure  time  series  did  not  reveal  a  sudden  jump  as  was  observed  in  pressure  sensors  1 , 
2,  and  3.  Its  cause  is  unknown  and  another  potential  source  of  error. 
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VI.  CONCLUSIONS 


The  data  analyzed  for  a  barred,  steep  beaeh  with  low-slope  incident  swell  waves 
did  not  compare  well  the  calibration  given  by  BJ85.  Shoaling  was  erratically  predicted  in 
the  model.  Reflection  is  likely  a  source  of  error,  which  is  not  accounted  for  in  the  model. 
In  order  to  achieve  more  conclusive  and  informative  results  from  these  data,  they  should 
be  filtered  so  only  those  at  high  tide  are  considered.  Additionally,  reflection  from  the 
steep  beach  should  be  considered  in  the  model. 
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